function y = euler(x0, a, n, y0, f)

	% Intrari:
		%x0=capatul stang al intervalului de integrare
		%a = lungimea intervalului
		%n = numarul de puncte
		%y0 = conditia initiala
		%f = functia de integrat y'=f(x,y)
	%Iesiri:
		% y = vectorul aproximatiilor solutiei
  	y = zeros( n+1, 1);
  	y( 1 ) = y0;
  	h = a / n;
  	for i = 2 : n+1
        	x = x0 + i * h;
        	y( i ) = y( i-1 ) + h * feval( f, x, y( i-1 ) );
        endfor
endfunction

